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CNj Abstract. This review highlights the properties of the planetary nebulae in external galaxies as 
^ tracers of hght, of the stellar population properties, and of the distances and kinematics of the 
^ \ parent galaxies. Recent results on the kinematics of the outer regions in giant elliptical galaxies 
. and on the luminosity specific PN numbers (the a parameter) in these systems are presented, 
based on current surveys of planetary nebulae with the Planetary Nebulae Spectrograph (PN.S) 
and other instruments. Finally a brief discussion is given of planetary nebulae as tracers of the 
diffuse light in the nearby clusters, such as Virgo and Hydra I. 
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Q^i 1. Introduction 

Q , Planetary Nebulae (PNs) are stars that evolve from the asymptotic giant branch (AGB) 
\-[ ' to their final destiny as white dwarfs. The majority of stars between 1 and 8 Mq evolve 
^ , through a PN phase. As a consequence, in old stellar populations as those observed in 
bulges and early-type galaxies, most stars will reach the post AGB phase and a fraction 
of those stars will go through a PN phase, before ending their lives as white dwarfs. 

In a PN, the central star emits most of the light in the UV part of the spectrum. 
The nebular shell of the PN is able to convert the UV ionizing photons into various line 
emissions from the UV to the optical and down to the NIR. Up to 15% of the UV emitted 
, energy by the central stars is re-emitted in the [OIII] A 5007 A line, the brightest optical 
' emission of a PN (Dopita et al. 1992). 

, We can integrate the whole [OIII] fiux emitted from a PN and derive the tosoo? mag- 
' nitude defined as 

m5oo7 = -2.51ogF5oo7- 13.74 (1.1) 

L| ' (Jacoby 1989). For a PN population observed in an external galaxy, we can derive the 
. ^ \ PN luminosity function (PNLF) for that population, which is empirically shown to follow 
' an analytical formula 

S ; ^(msoor) = C X e^-^o^™™" x [1 - e3(™*-"™''^)] (1.2) 

(CiarduUo et al. 1998) where to* is the apparent magnitude of the bright cut-off. This 
analytic formula combines the observed behavior at the bright end, which is believed 
to originate from the most massive surviving stellar cores (Ciardullo et al. 1989), and 
the slow PN fading rate caused by the envelope expansion at the faint end (Heinze & 
Westerlund 1963). When this equation is integrated down 8 mag from to*, it provides 
the total number NpN of PNs associated with the total luminosity of a given galaxy. The 
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magnitude range of 8 mag from to* allows us to account for the faintest known Galactic 
planetary nebula (Ciardullo et at 1989). 

In external galaxies, PN are most often observed using a narrow band filter centered 
on the redshifted wavelength of the brightest emission line, the PN [OIII] A 5007 A line. 
In galaxies at distances larger than 2 Mpc, the PN emissions are spatially unresolved. 

Planetary Nebulae trace light. The PN population in an external galaxy is expected to 
trace light because the luminosity-specific stellar death rate (B) should be independent of 
the precise state of the underlying stellar population (Renzini & Buzzoni 1986). Buzzoni 
et al. (2006) computed the specific evolutionary flux B for simple stellar populations using 
different Initial Mass Functions (IMFs) and different ages, and measured a variation of 
B by at most a factor of 2. We can write the total number of PNs, Npiy, associated with 
a given bolometric luminosity from a simple stellar population as 

NpN = B X Ltot X TpN (1.3) 

where 

• NpN is the total number of PNs associated with the parent stellar population. NpN 
is obtained by integration of the PNLF 8 mags down the bright cut off; 

• B is the luminosity-specific stellar death rate or specific evolutionary flux; 

• Ltot is the total bolometric luminosity of the parent simple stellar population; 

• TpN is the PN visibility lifetime, i.e. the time for the nebula to be detectable in 
either [OIII] or Ha surveys. The upper limit of tpn is 30000 yr, that is the time required 
for a PN nebular envelope to fade entirely because of its expansion (Heinze & Westerlund 
1963). 

We can define the luminosity specific PN number, a-parameter for short, as 

« = - B X rpjv (1.4) 

which provides us with a population averaged-measurement of the PN visibility lifetime 
Tp N for a PN population. Because of the small variations of the specific evolutionary fiux 
B, variations of the observed a values indicate primarily different values of Tpjv for PNs 
associated with the corresponding stellar populations. 

Extended PN samples were collected for early-type galaxies (ETG) as part of the 
Planetary Nebulae Spectrograph (PN.S, Douglas et al. 2002) ETG survey 
(littp://www. astro. rug.nl/~pns/). In ETGs, the PN number density profile follows the 
surface brightness of the light distribution, as shown for the PN.S galaxies (Coccato 
et al. 2009) and for NGC 4697 (Mendez et al. 2001), M 60 (Teodorescu et al. 2011) and 
NGC 1399 (McNeil et al. 2010). Work is currently ongoing for a sample of SOs galaxies 
(Cortesi et al. 2011; A. Cortesi, 2012 Ph.D thesis). 

Planetary Nebulae trace stellar populations. Buzzoni et al. (2006) investigated the time 
evolution of the luminosity specific PN number - the a parameter - using template galaxy 
models along the entire Hubble morphological sequence E-Sa-Sb-Sc-Sd-Im. The results 
indicate that the a parameter values predicted from these models depend on the mass 
loss relation during the evolution from the turn-off of the zero-age-main-sequence to the 
post-asymptotic giant branch phase; see Figure 7 in Buzzoni et al. (2006). When the a 
values are computed for the empirically calibrated Initial Final Mass Relation (IFMR) 
of Weidemann (2000), there results a flatter distribution from Im to E, see Figure [TJ 

Figure [T] compares the a model values with measurements from PN samples in ETGs 
and in galaxies from the Local Group. While model predictions and observed values 
are consistent for stellar populations in galaxy morphological types from Im to Sa, the 
observed a parameters in ETGs with red and old stellar populations are smaller than 
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Figure 1. A comprehensive overview of the luminosity-specific PN number in Local Group 
galaxies (open star symbols) and external ellipticals (solid and open triangle symbols). Su- 
perposed on the plot, there are template galaxy models. Galaxy evolution along the whole 
E-Sa-Sb-Sc-Sd-Im Hubble morphological sequence is tracked by models from 5 to 15 Gyr, with 
the latter limit marked by the big solid dots. Two model sequences are reported on the plot 
assuming an IFMR as from the standard case of a la Reimers mass-loss parameter ri= 0.3 (lower 
sequence), and from the empirical relation of Weidemann (2000) (upper sequence). An indica- 
tive estimate of the mean representative PN lifetime (in years) is sketched on the right scale, 
according to Equation 11.41 From Buzzoni et al. (2006). 
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Figure 2. Correlation between the oib,i.o parameter and FUV— V color, measured from total 
extinction-corrected magnitudes FUV (from GALEX) and V (from RC3), from Coccato et al. 
(2009). Filled circles: sample galaxies for which a is calculated from the PN.S ETG sample data. 
Open diamonds: a values taken from Buzzoni et al. (2006). 

the predicted values by up to a factor 7, see also Hui et al. (1993) and Ciardullo et al. 
(2005). 

The relative PN deficiency in ETG supports the presence of a fraction of low-mass 
cores with Mcore ^ 0.55Mq: for such low mass cores, rp^r may become shorter because 
the time required for the excitation of the nebular envelope increases, and by the time it 
happens, the nebular shell may be near to being lost. Furthermore the part of the stellar 
population with Mcore ~ 0.52Mo may omit the PN phase entirely. These evolved stars 
may provide an enhanced contribution to the hotter horizontal branch (HB) and post 
- HB evolution as directly observed in M 32 and in the bulge of M 31, which in turn 
implies that the far UV fiux in ETGs should anti-correlate with the measured a values. 
Such an anti-correlation is indeed observed in the PN.S ETG sample, see Figure [2j 
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A possible external mechanism that can act on a PN to shorten its rp^v is the inter- 
action with the hot intracluster medium. The ram pressure stripping by the hot high 
density gas in brightest cluster galaxies (BCGs) coupled with the large velocity of a PN 
through the medium can reduce the material in the PN envelope, therefore effectively 
reducing tpat (Dopita et al. 2000, Villaver & Stanghellini 2005). 

Planetary Nebulae populations trace distances. The luminosity function of the [OIII] 
magnitudes associated with a PN population can be described by the simple analytical 
formula given in Eq. 11.21 whose absolute magnitude of the bright cut-off M* is observed 
to be invariant in different morphological types, from late to early types (CiarduUo et al. 
2002). Because of the presence of a bright cut-off, the PNLF can be effectively used as a 
secondary distance indicator in both ETGs and late-type galaxies. One can ask whether 
the invariance of the PNLF bright cut-off for populations of different ages and metallicity 
can be explained on the basis of the evolution of simple stellar populations. 

From single star evolution, we know that a PN's peak flux in the [OIII] A 5007A line 
is proportional to its core mass (Vassiliadis & Wood 1994), the core mass is proportional 
to the turn-off mass (Kalirai et al. 2008), and the turn-off mass decreases with increasing 
age of the stellar population (Marigo et al. 2004) . Therefore the absolute magnitude of 
the bright cut-off M* should become fainter in older and evolved stellar populations. The 
observed invariance of the PNLF bright cut-off indicates that the brightest PN may not 
be the end result of simple evolution of single stars, but rather have different progenitors, 
including blue stragglers (CiarduUo et al. 2005) and symbiotic stars (Soker 2006). 

2. Planetary Nebulae as kinematical tracers 

In addition to measuring the TO5007 magnitude and determining the distance to a given 
galaxy via the PNLF, the [OIII] A 5007 A emission allows us to measure of the line-of- 
sight (LOS) velocity of a PN. By measuring vlos for the PN population in the outer 
regions of galaxies we can obtain the kinematics of the stellar population at large radii 
and derive the projected two-dimensional velocity field. 

Since the work of Arnaboldi et al. (1994) on NGC 1399 and of Hui et al. (1995) on Cen- 
taurus A (NGC 5128), there has been a considerable success in mapping the kinematics 
of the stellar population in the outer regions of ETGs out to 5i?e (Arnaboldi et al. 1996, 
Arnaboldi et al. 1998, Mendez et al. 2001, Romanowsky et al. 2003, Coccato et al. 2009, 
Napolitano et al. 2009, McNeil et al. 2010, Napolitano et al. 2011, Teodorescu et al. 2011, 
Cortesi et al. 2011), and intracluster light (ICL) in nearby clusters (Arnaboldi et al. 2004, 
Gerhard et al. 2005, Gerhard et al. 2007, Doherty et al. 2009, Ventimiglia et al. 2011). 

In the galaxies studied thus far, the PN density and kinematics is everywhere consistent 
with the integrated light measurements within the errors (Coccato et al. 2009), with the 
possible exception of NGC 4697 (Sambhus et al. 2006). The outer halos kinematics from 
the PN samples in ETGs show the following properties: 

• there is a dichotomy in the outer halo kinematics, the velocity dispersion profile is 
either slowly falling or rapidly falling (Coccato et al. 2009) with increasing radius. 

• The rotational properties of the halos traced by the ratio v/a where v is the circular 
velocity and a is the velocity dispersion correlates with that within i?e for the greater 
part of the sample galaxies observed so far. 

• The slow/fast rotators division in the outer halos (Coccato et al. 2009) is similar as in 
the cores (Emsellem et al. 2007). There are some more complicated cases, as kinematic 
misalignment is more frequent at large radii (Shih & Mendez 2010, Teodorescu et al. 
2011). 

• A fraction of morphologically-undisturbed ellipticals show large distortions of their 
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kinematics or presence of kinematic sub-components at large radii, as in the outer halos 
of NGC 1399 (McNeil et al. 2010) and NGC 3311 (Ventimiglia et al. 2011). 



3. The brightest Planetary Nebulae in the Virgo Cluster 

The PN populations in the core of the Virgo cluster are studied via narrow band imag- 
ing surveys, using the "on-band, off-band" technique, in which a field is imaged through 
a narrow band filter, centered at the redshifted [OIII] A 5007 A emission in these clusters, 
and a broad band filter. PN candidates are those unresolved sources that have a color 
excess in [OIII] and no detected continuum (Arnaboldi et al. 2002). Tens of intracluster 
PN (IPN) were discovered in the Virgo cluster, in its core and in outer targeted fields 
(Arnaboldi et al. 2002, Okamura et al. 2002, Arnaboldi et al. 2003, Feldmeier et al. 2003, 
Feldmeier et al. 2004, Aguerri et al. 2005, Castro-Rodriguez et al. 2009). Such surveys 
are complete 0.5 mag down the PNLF, therefore they sample mostly the brightest PNs 
in these stellar populations. The spectroscopic follow-up with multi-object (Arnaboldi 
et al. 2003) and multi fiber spectrographs (Arnaboldi et al. 2004, Doherty et al. 2009) 
allow 

• to detect the second [OIII] A 4959 A emission besides A 5007 A in single spectra, 

• measure the radial velocities and associate PNs with the dynamical components 
along the LOS, 

• determine PN physical parameters, the rnsoo? and the velocity vexp derived from 
the FWHM of the [OIII] A 5007A emission, i.e. the expansion velocity of a PN nebular 
envelope. These PN physical parameters can be correlated with those of its parent stellar 
population (Arnaboldi et al. 2008). 

In case of M87, the BCG in the Virgo cluster, the PN LOS velocity distribution traces 
the M87 halo out to 150 kpc and shows that the halo PN population has a slightly brighter 
cut-off than that in the galaxy's central region, by 0.3 — 0.4 mag, see the observed M87 
PNLF in Figure |3l Because the LOS distance to these PNs is known, via the physical 
association with the M87 halo, their absolute magnitude is also known, and so are then 
the core masses required to power the detected [OIII] A 5007 A fluxes: for the brightest 
PN in the M87 halo, the core masses are ~ O.62M0. The IFMR for solar metalhcity 
stars predicts that PN progenitors with 2.2AfQ give final core masses of 0.62AfQ. Turnoff 
masses of ~ 2Mq belong to populations with ages ~ IGyr. It is unlikely that such 
progenitors exist in the Virgo core, as the stellar population is mostly old (Williams 
et al. 2007): therefore close binaries and blue stragglers are indicated as most likely 
progenitors for these intrinsically bright PNs. 



4. The brightest Planetary Nebulae in the Hydra I Cluster 

The Multi-Slit Imaging Spectroscopy Technique. The nearest clusters beyond Virgo and 
Fornax are at distances of about 50 Mpc, and the brightest PN in the Hydra I and the 
Centaurus clusters have fiuxes of 8 x 10""'^^ erg s~^ cm~^. They cannot be detected using 
standard narrow band imaging because the noise in the 40 — 60 A sky is of the same 
order of the signal we want to detect. 

A successful technique to detect the [OIII] A 5007A emission from the brightest PNs in 
these clusters is the "Multi-Slit Imaging Spectroscopy technique" (MSIS, Gerhard et al. 
2005, Arnaboldi et al. 2007). This is a blind search technique that combines a mask of 
parallel multiple slits with a narrow band filter, centered at the redshifted [OIII] emission 
of the PNs in these clusters, plus a dispersing element to obtain spectra of all PNs that 
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Figure 3. Left panel: the observed PNLF in the Virgo cluster core fields, in 0.3 mag bins. 
Long dashed blue line: PNLF for the entire spectroscopically confirmed sample of PNs in the 
Virgo cluster core. Red continuous line: PNLF for the M87 halo PNs. Short dashed black line: 
analytical formula for the PNLF with m* — 26.0. Right panel: LOS velocity distribution for the 
14 M87 halo PNs. The systemic velocity of M87 is also indicated. From Arnaboldi et al. (2008). 



lie behind the slits. The sky noise at the emission line of a PN conies from a few A only, 
depending on slit width and spectral resolution (Gerhard et al. 2005). 

PN population in the halo of NGC 3311. The MSIS observations of the Hydra I cluster 
carried out with the F0RS2 spectrograph at the VLT detected 56 PNs and 26 background 
galaxies in a 6'. 8 x 6'. 8 field centered on the BCG NGC 3311 in the cluster core (Ven- 
timiglia et al. 2011). In addition to the study of the halo kinematics, Ventimiglia et al. 
(2011) measured the luminosity specific PN number a for the stellar population in the 
extended NGC 3311 halo, and obtain a value that is a factor 4—6 lower than the predicted 
value of log a — —7.30 from the a vs. FUV — V color correlation for the ETC populations 
shown in Figure [51 In Figure 0] the observed cumulative PN number for the NGC 3311 
halo (green line) and the total NGC 3311 plus ICL light (black line) in the Hydra I core 
are plotted together with the predicted cumulative PN number for the NGC 3311 plus 
ICL (using logo; = —7.30, red line), as function of the radial distance from the center of 
NGC 3311. The large discrepancy between the red and green lines indicates a decrement 
of PNs in the NGC 3311 halo. This PN deficit may signal that additional effects are at 
work that shape the properties of the PN population in the NGC 3311 halo. 

A possibility is that ram pressure against the hot X-ray emitting gas in the halo 
of NGC 3311 is high enough to severely shorten tpn, the visibility time of a PN, see 
Section [TJ With a density of the intracluster medium (ICM) inside 5' around NGC 3311 
of about picM — 6 X 10^"^ cm~^ and a typical velocity oi v — -\/3 x 450 kms^^ ^ 800 
kms^^, the ram pressure on the NGC 3311 halo PNs is ^ 40 times stronger than in the 
simulated case by Villaver & Stanghellini (2005) for the Virgo cluster. The ram pressure 
effect could be much stronger in the Hydra I cluster and shorten significantly rpjv for 
the NGC 3311 halo population. 
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Figure 4. Observed and predicted cumulative PN number for the NGC 3311 halo and the total 
(NGC 3311 plus ICL) light in the Hydra I core, as function of the radial distance from the center 
of NGC 3311. Red line: predicted cumulative number of PNs computed using logo = —7.30. 
Green line: measured cumulative number of PNs associated with the NGC 3311 halo on the 
basis of the vlos- Black line: measured cumulative number of PNs corresponding to the total 
light in the Hydra I core. From Ventimiglia et al. (2011). 

5. Summary and Conclusions 

In this review we illustrate that PN populations are ubiquitous to the light associated 
with external galaxies and the diffuse light in the cores of nearby clusters. Deep photom- 
etry and kinematics of PN populations show that they trace star light and motions in 
the extended, luminous halos of BCGs and ICL. 

The PN population shows dependencies on the stellar population age and metallicity. 
Recent observations of PNs in the Hydra I cluster indicate that there may be an additional 
effect by the hot ICM, which could decrease the visibility lifetime tpn of PNs. Further 
theoretical and observational work is needed to quantify effects of stellar population age, 
metallicity and ICM on both the PNLF bright cut-off and the luminosity specific PN 
number a. 

The study of the PNs in the ICL and the halos of galaxies is important for constraining 
how and when these components formed. In the case of the Hydra I cluster the results 
indicate that the build-up of the BCG halo and the nearby ICL continues, and that the 
formation of the outer halo and ICL is an on-going and long lasting process. 
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